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A special class of 5 = 1 spin ladder hamiltonians, with second- neighbor exchange in- 
teractions and with anisotropies in the z-direction, can be mapped onto one-dimensional 
composite S — 2 (tetrahedral 5=1) models. We calculate the high temperature expansion 
of the Helmoltz free energy for the latter class of models, and show that their magnetization 
behaves closely to that of standard XXZ models with a suitable effective spin S e ff, such that 
5 e //(l + S e ff) — (Sj), where Si refers to the components of spin in the composite model. 
It is also shown that the specific heat per site of the composite model, on the other hand, 
can be very different from that of the effective spin model, depending on the parameters of 
the hamiltonian. 
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Many materials have been described by low-dimensional spin models 1] such as XXZ, ladder, 
tetrahedral, dimmer chain, and mixed spin models. Different predictions are obtained from 
one-dimensional and quasi-one-dimensional models that can be verified experimentally; e.g., 
the 5 = 1/2 antiferromagnetic Heisenberg chain is a gapless model, whereas the even-legged 
antiferromagnetic Heisenberg ladder model has a gap in its energy spectrumQ]. Attention has 
also been drawn to composite spin models; for instance, the T = phase diagram for the 
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5 = (l/2©l/2) model has been studied by Solyom and Timonen This composite spin model 
is equivalent to the tetrahedral 5 = 1/2 model - applied to the study of the properties of the 
tellurate materials C^T^OsCZs and Cu2Te20^Br^. On the basis of experimental results it was 
argued that these materials could be appropriately described by the noninteracting tetrahedral 
5 = 1/2 modelfl. In Ref. |a| we obtained the high temperature expansion of the Helmholtz 
free energy of the tetrahedral 5 = 1/2 model up to order /3 5 . 

A common feature shared by all composite spin models is that the modulus of the spin at 
each site of the chain is not constant; instead, one has a random distribution throughout the 
chain. The presence of spin 5 = in a given site of the chain can be interpreted as nonmagnetic 
impurity at that sitejfj]. 

By the appropriate design of molecules, it is possible to obtain a variety of spin systems. One 
example is the organic compound BIP-TENOpj, found to be a 5 = 1 spin ladder. At T = 
this compound exhibits a plateau- like anomaly at 1/4 of the saturation magnetization, which 
has attracted much attention due to its full quantum nature. In Ref. i| a quasi one-dimensional 
spin model with second- and third-neighbor antiferromagnetic exchange interactions suitably 
explains the mechanism of the 1/4 magnetization plateau of this material; it does not include 
any anisotropy in the z-direction, though. Numerical analysis has been applied to study its 
phase diagram 9] at T = 0; even in the absence of third-neighbor interactions the magnetization 
plateau is shown to be present. 

In Ref. ^(| we presented a method to calculate the coefficients of the cummulant expansion 
of the Helmholtz free energy, in the thermodynamic limit, of any chain model whose hamiltonian 
satisfy periodic boundary conditions, possess invariance under spatial discrete translations and 
include interactions between nearest neighbors. 



By a suitable choice of constants in the hamiltonian of Ref. 8|, the 5=1 ladder with second- 

n 

neighbor interactions (5 = 1 tetrahedral model, cf. model B in Ref.jllJ; at each site, spin- 1/2 
degrees of freedom are replaced by those of a fundamental spin- 1) can be mapped onto a family of 
one-dimensional composite 5 = 2 chain models. In the present work, we increase the parameter 
space of the model in Ref.|8| by introducing anisotropies in the z-direction, calculate the high- 
temperature expansion of its Helmholtz free energy and compare its thermodynamics, in the 
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same regime of temperature, to the standard (irreducible representation ) spin- S XXZ model. 
The high-temperature expansion (or /3-expansion, where = 1/kT, k is the Boltzmann constant 
and T is the absolute temperature) of the Helmholtz free e nerg y for the composite S = 2 XXZ 
model is obtained by applying the method presented in Ref. 10|. From the analytical results we 
may obtain the thermodynamic functions of the model at high temperatures; in particular, the 
magnetization M as a function of the external magnetic field h and the specific heat per site Cl 
as a function of (3. 

The hamiltonian of the quasi-one-dimensional spin model is 



N ! 

i=l 

+ J[(cTj, cr i+ i) A <8> 1 T + (o"i, T i+ i) A + (n, a i+ i) A + 1 CT <8> (r i5 r i+ i) A ] 
- h{a\ ® 1 T + l a ® rf)}, 

and it is subject to periodic boundary conditions. We use the same notation as in Ref. 
{Ai,B k ) A = A x l ®B x k +A y l ®Bl + AAf B z k , with A t = {A X ,A\, Af) and B k = (B x k , B y k , B k ), 
introducing the anisotropy in the z-direction. For Ao = 1 and A = 1, ffl equals to the sum of 

n 

hamiltonians (2) and (3) of reference |8| for the S = 1 ladder with second-neighbor exchanges 
(J3 = 0), for the special case J\ = J2- The distinct 5 = 1 variables <jj and t% are related to the 
p- and r-lines of the dumb-bell, respectively (cf. Fig.l of Ref.[y]). 

We define the composite spin operators Si at the i-th site as Si = Oi (8) 1 T + l a (8) 7^. Here, 
1 CT and 1 T are the identity operators in a- and r-space, respectively. Using the composite spin 
Si, the tetrahedral S = 1 hamiltonian is rewritten as a composite S = 2 chain hamiltonian 

Kq-id = £[- 2J ol + 0(Si,Si)i + j (S+Sr +1 + SrS+ 1 + ASfSf +1 ) 

i=i 

-hSi + d(Sff], (2) 

where 5 = d = ^(Ao — 1), Sf = ^{Sf i ^"Sf ), and 1 is the identity operator, represented 
by a 9 x 9 identity matrix. The block matrix representations of the composite spin operators in 
©, in the basis of eigenstates of and Sf, are 
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We point out that the operator is a constant of the motion, although the hamiltonian 
can be interpreted as mixture of 3 kinds of spin {S = 0, 1 and 2), randomly distributed along 
the chain, with its probability depending on the constants of the hamiltonian and on the tem- 
perature. 

An interesting limit of the hamiltonian (J2J) is g = Jo/2 — > 0, for finite values of d = Jo ^° . In 
this limit it acquires a single-ion anisotropy term. 

As in Ref.ja|, the presence of a spin S = in a site of the chain can be interpreted as 
nonmagnetic impurity at that site. The model thus encompasses the presence of randomly 
distributed impurities along the chain, without hindering the application of the method of Ref. 
[lol |. since the hamiltonian © has only nearest-neighbor interactions and is invariant under 
spatial translations. 

Appendix presents the /3-expansion of the Helmholtz free energy Wq-w of the composite 
S = 2 XXZ model, up to order /3 6 . Although this expansion has a large number of terms, it is 
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easily differentiaded by any symbolic computer language, yielding the thermodynamic functions 
of the tetrahedral 5 = 1 model. 

For finite temperatures the results are insensitive to the sign of J. If we redefine the param- 
eters of the hamiltonian © in units of J, the expression (|Alj) of Wq—xd becomes an expansion 
in {(3 J). For the sake of simplicity, we choose J = 1 in what follows. 

As we are discussing the case of the quasi one-dimensional model (the tetrahedral S = 1 
model) onto a chain model, it is interesting to compare the high-temperature thermodynamic 
properties of the latter, to those of the standard spin-S" XXZ model, for distinct values of S. 

The free energy Wq-w and the mean value of the squared norm of spin (S 2 ), for the S = 1 
ladder model with second- neighbor exchange interactions, are related by (S?) = for 
g = 4r and d = 4r(Ao — 1). (S 2 ) is a function of temperature and of the parameters of the 
hamiltonian ©. At = 0, the universal value (S 2 ) = 4 is obtained, since in this limit of 
temperature we have N independent spins with equal probability to be at S = 0, 1 or 2. 

FigE shows (S 2 ) as a function of j3 for two sets of values for the parameters of the composite 
model. (S?) is not very sensitive to the value of h, within the range of h where its /3-expansion 
is a good approximation. We let h = in Fig^ m Fig^, we have A = —0.3,g = 0.5 and 
d = —0.35; and in Fig^3 A = l,g = —0.5 and d = 0. For the sake of comparison, Fig^also 
shows straight horizontal lines that correspond to the squared norms of spin per site for the 
standard XXZ model with S = 3/2 and 2. 

FigEl shows the entropy S as a function of /3, in units of Boltzmann's constant (S = Z? 2 ^-), 
for the same set of parameter values as in Fig^ hi both cases the entropy of the tetrahedral 
S = 1 model is higher than the entropy of the composite S = 2 XXZ model. This result agrees 
with the fact that the number of degrees of freedom of the composite S = 2 model is higher 
than that of the standard (irreducible) S = 2 XXZ model. 

In reference Rojas et al. derived the high temperature expansion of the Helmholtz free 
energy Wg of the standard spin-S XXZ model, for arbitrary value of S (S being semi-integer 
or integer), up to order /3 6 . From their result, we obtain the expansion of the magnetization 
(M) of this family of models (M = — ^j^)- Appendix iBl shows this expansion up to order /? 6 . 
Postulating the validity of expansion (|B1|) for real positive values of S, we define an effective spin 
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FIG. 1: Expectation values of the squared norm of spin (Sf ) as a function of (3 = 1/kT for the composite 
5 = 2 model (solid lines) and the standard spin-,!? model, with 5 = 3/2 (dashed lines) and S = 2 (dot- 
dashed lines), for a vanishing magnetic field h. The values of parameters of the hamiltonian @ have 
been chosen as (a) A = —0.3, g = 0.5 and d = —0.35 and (b) A = 1, g = —0.5 and d = (no anisotropy 
in the z-direction). 

value S e fj such that 5 e //(5 e // + 1) = (S^) where (S?) relates to the composite 5 = 2 model. 

Obviously, the effective spin depends on the values of the parameters of the theory (see Fig^) 
and, in general, it is neither integer nor semi-integer. FigEk compares the magnetization of the 
tetrahedral 5 = 1 model and that of the effective XXZ model at (3 = 0.2, for A = l,g = —0.5 
and d = 0. FigEJa shows the percent error of the two curves in Fig|3K- Such error is less than 
1.3% for h e [0,1.4]. 

This similarity of behavior is not shared, in general, by the specific heat per site Cl {Cl = 
—f3 2 d ffl ) , as shown in Fig@J 
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FIG. 2: Comparison of the entropy 8 as a function of j3 = 1/kT for the composite S = 2 XXZ model 
(solid line) and the irreducible S — 2 XXZ model (dashed line). The values of parameters in Figs [2^ and 
03 are the same as in Figs^ and^ 3 , respectively. 

Fukushima et al\\\ calculated the /3-expansion of the specific heat per site of the spin-S 1 XXZ 
model, for arbitrary (semi-integer or integer) values of S, up to order f3 . Proceeding similarly 
as done for the magnetization, we postulate the validity of their expansion of the specific heat 
per site to any real positive value of S. Fig@^ shows the specific heat per site as a function of 
(3 for the tetrahedral S = 1, the standard S = 3/2 and the effective S e ff XXZ models; parameter 
values h = 0, A = —0.3, g = 0.5 and d = —0.35 are the same as in Fig^,. At high temperatures 
the curves are close: for (3 G [0.13,0.27] the percent error is smaller than 5%, but for (3 ~ it 
goes up to 27%. 

FigEb was obtained for the same values of parameters used in Fig^^ that is, h = 0, A = 1, 
g = —0.5 and d = 0. The Cl curves for the effective and the composite 5 = 2 models are 



8 



p = 0.2, A = 1, g = -0.5 and d = 




FIG. 3: The solid line in FigEt shows the magnetization M as a function of the external magnetic field 
h at p = 0.2 and for the following values of parameters of the composite S = 2 XXZ model©: A = 1, 
g = —0.5 and d = 0. The dashed line is the magnetization curve for the effective S e ff XXZ model. Figl^ 
shows the percent error of the two curves in Fig|3Ji as a function of h. 

very different: the percent error for ft € [0, 0.2] varies from 20% at (3 = 0, to 41.8% at /3 = 0.2. 
Although there are 5 = and S = 1 sites in a composite 5 = 2 XXZ chain, its specific heat is 
larger than that of a standard 5 = 2 chain, for /? > 0.15. 

In summary, we applied the method of Ref.^^] to calculate the high temperature expan- 
sion (/3-expansion) , up to order /3 6 , of the Helmholtz free energy of a special class of quasi 
one-dimensional models, the 5 = 1 ladders with second neighbor exchange interactions (the 



tetrahedral 5 = 1 models 



We have increased the parameter space of this class of hamil- 



tonians by introducing anisotropies along the z-direction. This class of quasi one-dimensional 
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FIG. 4: (a) The specific heat per site Cl as a function of j3 for the composite S = 2 model (solid line), 
the S = 3/2 XXZ model (dot-dashed line), and the effective S e ff XXZ model (dashed line), for h — 0, 
d = —0.35, A = —0.3 and g = 0.5. (b) The same, for the composite S = 2 model (solid line), the standard 
S = 2 XXZ model (dot-dashed line) and the corresponding effective irreducible spin S e f / model (dashed 
line) for h = 0, d = 0, A = 1 and g = —0.5. (The values of S e ff are different, in each case.) 

models can be mapped onto one-dimensional composite S = 2 XXZ models, with spins S = 0, 1 
and 2 randomly distributed along the chain. At high temperatures, the thermodynamic prop- 
erties of the composite S = 2 XXZ model and the standard spin-5 XXZ model are expected to 
differ, in principle; this is the case of the specific heat. However, a somewhat surprising result 
is that the magnetization of the tetrahedral 5 = 1 model follows that of an effective irreducible 
spin S e ff model very closely (up to a percent error smaller than 2%). The effective spin is such 
that S e f(S e f + 1) = (S?), and (Sf) is the expectation value of the squared norm of spin per site 
of the composite S = 2 XXZ model. The function (Sf) is a real- valued continuous function; 
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S e ff may not be either integer or semi-integer. 
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APPENDIX A: THE HELMHOLTZ FREE ENERGY OF THE TETRAHEDRAL S = 1 

MODEL 

The following expression for the high-temperature expansion up to order /3 6 of the Helmholtz 
free energy refers to the composite 5 = 2 model (that is, the tetrahedral 5 = 1 model), whose 
dynamics is driven by the hamiltonian ©. It is valid for any real values of Jo,g, J, A, d and h. 
(Observe that these values may be positive, null or neg ative). When g = 4f and d = 4?(A - 1), 
it yields the free energy of the hamiltonian Q , for the 5 = 1 ladder model with second- neighbor 
interactions. 



Wq-id = -^-p + Wo + W 1 (5 + W 2 P 2 + W 3 /3 3 + W 4 /3 4 + W 5 /3 5 + W 6 (3 6 + 0(/f), (Al) 
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APPENDIX B: HIGH TEMPERATURE EXPANSION OF THE MAGNETIZATION OF 

THE SPIN-5 XXZ MODEL 

From the high temperature expansion of the Helmholtz free energy Ws of the irreducible 

n 

spin-S" XXZ model in Ref. 12], we obtain the magnetization M as function of the squared norm 
of the spin s(s + 1), up to order (3 6 : 



M = = M 1 (3 + M 2 f3 2 + M 3 /3 3 + M 4 (3 4 + M 5 (3 5 + M 6 /5 6 + 0(/3 7 ) (Bl) 
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Mi = k -^±^ (B2) 



/ 4 , , 2 . , \ 2 / ,,2 hds (s + 1) ,-,„. 
M 2 = - ( — hd+-Ah\ s 2 (s+l) 2 + y - >- (B3) 



M3 = -(l| /l -9i5 /ld2 -li A2/l -l| A/id ) s3 ( S + 1 ) 3 

— h+— hd 2 + — h 3 + — A 2 h+ — Ahd) s 2 (s + l) 2 
30 105 45 45 45 

- f-^^ + go^) s(s + 1) (B4) 



^--(-44 fed+ S A ^ + ^ A2 ^ + 2^ A3 ^4^-l^^) S4 ^ + 1 ) 4 

_ (-A. Ah 3 - i51 Ahd 2 - — ft 3 d- — hd 3 - — A 3 ft 
V 135 4725 945 4725 675 

-S^^-l^^-^^)^^^^^^^^^^^ 

+ 94 ft3d+ 7^ A ^-l4 A ^ + lf5 A2ftd+ lf5 A/ld2 ) s2 ( S + 1 ) 2 

-(l4 ft3d -90^ 3 ) s(s + 1) (B5) 



32 . , , 64 . 3 , , 718 . 4 , 
■ A h d - — — A hd~ — — — A h 



v 2835 1701 42525 

4 h-JLA hd >-3,A>h d > + ^A> h+ JLh^-JLhd 2 



2835 42525 2835 42525 93555 14175 

s 5 (s + 1) 5 - ( *—hd 2 -^-hct l + — h 3 d 2 - — h-— h 3 

y ' \ 14175 31185 1575 567 525 

128 ., , 176 , s , , 352 , 3 64 . 2 , j2 32 . , 3 , 

H Ahd-\ A hd-\ A h d H A hd -\ Ah d 

14175 2835 14175 945 525 

383 a 2 , , 229 a4i 358 , 2 , 3 \ 4 , ,, 4 / 11 , , 2 184 , 4 



+ 42525 A fe+ 14175 A ^ + 4725 A h ) S > +1 > -(-1575 *" " 31185 hrf 

14175 210 14175 945 1575 567 

652 ., , 3 8 . 2 , , 2 136 . , 3 , 271 . 2 , 64 . 4 , 
A /id - — — A 2 hd - , Ah d— — — — A 2 h - ——— A h 



14175 135 14175 14175 14175 

+ J?4. A 2 h A g 3 ( 1)3 _ ( J}_ hd 2 1_ 4 _ _U_ 3 2 J_ _ J7_ 3 

14175 / v ; \6300 77 3150 840 6300 

-io h5 + ^ Ahd+ i- 5 Aahd+ w9 Ahda + h A2hd2 ~§k Ah3d 

+ is A2h+ lL A4h -^ A2h3 ) s2 ^ 2 

-W2 hd *-^ + m hid2 )^^ 
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( 1592 7 7 

Me = -(-1488375 ^ 

+ S h ^ + 638512875 * ^ " ^ h * + sS A ° * ^ 

35072 AI l4 70424 a4 , , 2528 . , 1672 a3i 5132 a5i 
■ A ft <T + -—— A 4 ft d + —— Ah- ——— A h + — —— A ft 



49116375 4465125 893025 178605 893025 

5792 a2uSS 1856 ^ 6 6 / ^ 304 3 

+ 1964655 A ^ + 1488375 Aftd J S (s + 1} " (,"1488375 hd ~ 606375 ftd 
+ JL fe 3 d+ 24512 ftd .__g^ ftV _^g_ Aa jggg, 
4725 70945875 66825 893025 33075 

44416 . , , 4 5942 . 4 , , 7628 . , 2 . 3 , 64 . , 3 
A 4 ft d + Ah+ — — A ft + — — A h J 



16372125 165375 893025 6615 2835 

2656. A5ft _ 9^ A 3 ft3 _ 128 A ,3 d2 _ 1328 A2ft 3 d _ ^ 



297675 14175 4725 14175 654885 

165375 J v ' V 992250 606375 675 

13136 5 J0gL ft 3 d , + _L h B d+ igg7_ Aahd+ jg4g. A 3 htfa 
70945875 467775 1575 59535 33075 

112864 AI , 4 13882 a4 , , 1543 17137 a3i 139 . , 3 

H A ft d H A ftd H A ft H A ft H A h 

5457375 496125 1190700 1190700 4725 

+ ^A 5 ft--^-A 3 ft 3 + -|-Aft 5 + -|l-A/ l a d 2 + -^A 2 / l 3 d 
99225 42525 175 1575 525 

^A^hd^^AhdA s\s + lf-(-^hd-^^h d * 
13475 165375 ) \ 220500 1819125 

+ 4^' l3d -7li^^ 5 -^ ft3d3 + 4^ ft5d -3i| A2/ld 

2606. A 3 fc ^_^ Ah ^_ 8111 A4fcd _ Aft _ «1 A8ft 



99225 5457375 992250 33075 66150 

+ Jl 2 _ A/i 3_^ A5 , + J! 2 _ A 3,3 + JL A , 5 + ^L Aft 3 d2 

14175 1323 14175 4725 14175 

+^A 2 ft 3 d- -^-A^h^-^AhdAs^s + lf-f^hd 
14175 1091475 110250 J y ' \ 88200 

+ -gg- Ad »-JL ft 3 d+ 177571 h1 ?-J*- h :t>-J»- h ' d 

161700 6300 18918900 6237 18900 

4 a 2 , , 74 , 3l , 2 4553 A , , 4 37 A 4 , , 13 A , 

H A h d H A ft d H A ft d H A ft d H A ft 

1225 33075 363825 66150 14700 

17 q 11 , 1 2 q q 11 4 q , 

+ 22050 A ^ 3150 A?1+ 26460 A ^ 2835 A >h + 4725 A ^ " 189 A ^ d 
-^ A2ft3d+ Jl5 A2 ^ 3+ l470 A ^ 2 ) s2(S + 1)2 
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